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ABSTRACT 

This article gives an brief outline of the applications of generalized special 
functions such as generalized hypergeometric functions, G-functions and H- 
functions into the general area of nuclear energy generation and reaction rate 
theory such as the energy generation in a simple stellar model and nuclear 
reaction rates in non-resonant and resonant as well as screened non-resonant 
situations. 

1 Introduction 

During the last few years several authors have been trying to get explicit and 
exact analytic representations for various problems in astropyhysics. Ana- 
lytical solutions are often extremely difficult to obtain and often lead to 
complicated intractable mathematical expressions. Due to these facts some 
researchers go for approximations, numerical results and computer print- 
outs. But such numerical results usually do not give much insight into what 
is happening between the stage of the formulation of the problem and the fi- 
nal computer printout. When the experimental results do not agree with the 
computer printouts resulting from a proposed theory one needs explicit ana- 
lytic formulation to modify the theory. Hence there has been renewed interest 
in getting explicit and exact theoretical results rather than numerical results. 
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During the past few years the authors have demonstrated that the techniques 
of generahzed special functions combined with statistical techniques are very- 
powerful in getting analytical results for many problems which were believed 
to be intractable mathematically. This paper will confine itself to the appli- 
cation of such techniques for studying energy generation and thermonuclear 
reaction rates for stellar models. 



2 Energy Generation Rate in Simple Stellar 
Models 



Some simple stellar models for looking into the internal structure of stars have 
been given by Kourganoff (1980) and Haubold and Mathai (1984, 1986b). 
Consider a spherically symmetric star in quasi-static equilibrium, that is, we 
exclude magnetic fields, rotation and other effects which are changing rapidly, 
so that we can take the star as a gas sphere in hydrostatic equilibrium. Let 
r denote the distance from the center to an interior point and R the radius 
of the star under consideration. 

2.1 Linear model 

The simplest model one can think of for the radial density p(r) is the follow- 
ing, 



p{r) = pc(l -s)- Pc{l -x),x^ r/R, p(0) = pe, p{R) = 0. (1) 



r 



In this case the mass 



M{t) = 47r / ep{t)dt = —Pcr\l 
Jo 6 



3x 

T 



(2) 



the pressure 




(3) 
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where G is the gravitational constant, and the temperature from the equation 
of state P(r) = kp{r)T{r)/MuH, where fi is the mean molecular weight, M„ 
is the atomar mass unit, k is the Boltzmann constant, is then 



36k 1 — X 

with Co = ci,C2 = — 19/5,C3 = 9/5, C3 = 9/5, c„ = for n > 4. The nuclear 
energy generation rate e(r) will depend on T(r) and p(r). If we take a simple 
model, 

e(r) = eo(po,To)(p(r)/po)"mr)/ro)^a > 0,/3 > 

then 

00 

eir)^d,il-xnY.CnXn^, (5) 

n=0 

where di is a constant. The luminosity L of the star is then, 

rR 

L — 4:7: r'^ p(r)e(r)dr 
Jo 

= d2 / X^il-xr+\Y.CnX''fdx, (6) 
•^0 n=0 

where ^2 is a constant. We consider a general integral 

gia,(3,j)= x^(l-xr+i(^c„x")^dx. (7) 



The L in (6) is evaluated in Haubold and Mathai (1984) for various structures 
of Z^^o ~ One such structure is when 

00 

^ Cnx"- = (1 + aix){l + a2x) . . . (1 + akx) 

n=0 

for a fixed k. In this case one can show that g{a,(3,j) can be evaluated in 
terms of a Lauricella function Fd- Then a particular case will go in terms of 
Gauss' hypergeometric function. 
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2.2 A general stellar model 

A linear decrease of the density from the center to the surface may not be 
very appropriate. Faster or slower decrease can be incorporated by using one 
more parameter and taking the model for the density as, 

Pir) = Pc(l-(^)') 

= p,{l-x'),x^r/R,p{0)^p,,p{R)^0. (8) 

In this case the mass, pressure, and luminosity will work out to be the fol- 
lowing. 

M(r) = ^ = p,i?V(l-^a;^). (9) 

, 1 S + 6 



(10) 



2 (5 + 2)(5 + 3) 2(5 + l)(5 + 3)' 



and 



L{R) = 4wo(Pc,7;)^ /^(3/^)-^(l-x)"+^+^ 

6 Jo 



2V. 1 (i + 2)(i + 3) + (i + l)(i + 3)/' 

For /3 = 1 or a positive integer the evaluation of the integral in (11) is not 
difficult. For the general case one can look at various situations as in Section 
2.1 and solve many of these in terms of generalized functions. 

2.3 A further generalized stellar model 

A model with more flexibility for the density is a two parameter model 

pW=Pc[i-(^)'r,5>o,7>o. 

In this case the mass will lead to the form 

t 



Mir) = 4npc Tt^ 
Jo 



^ , t r 

^^>0<-<-<l 



= r^-S^(-7,|^i;(^)1, (12) 



where 2-^1 is a Gauss' hypergeometric function. Pressure and temperature in 
this case can be obtained in series of hypergeometric functions. For example, 



P(r) 



^^Gpl \ ^ (-7)^(1 + 1) 



\ i'o^!(| + l + 7Uf + m)(f + m) 




- {-l)m[{^)\ 



(13) 



m!(f +m)(| + m) 



2F1 



(-74 + m;?+m+li(^r)}. 



(14) 



Note that when 7 is a positive integer then P{r) is available as a finite sum. 

3 Collision Probabilities 

In Haubold and Mathai (1984a) a detailed derivation of the collision prob- 
ability for nuclear reactions is given under the assumption that the nuclear 
velocity distribution remains Maxwell- Boltzmannian. The relation between 
the macroscopic nuclear reaction probability and the microscopic nuclear re- 
action cross section on the basis of equilibrium-thcrmodynamic arguments 
is studied in this paper as well as in Haubold and Mathai (1985). A ba- 
sic integral to be evaluated to compute the reaction probability, under this 
approach, is the following: 



More complicated integrals will appear when dealing with thermonuclear 
reaction rates in various non-resonant and resonant cases, see for example, 
Haubold and Mathai (1986, 1986a,b,c). First we will introduce a statistical 
technique of tackhng integrals of the type (15). Let us consider a shghtly 
more general integral 




(15) 



Jo 



, — n/m 



dy. 



(16) 



Theorem 3.1. 
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= p-'-(27r)(2-"-'")/2m^/2n(i-2'^'-)/2 (17) 

for i?(p) > 0, R{z) > 0, G{.) is a G-functfon and R(.) denotes the real part of 
(.). For a discussion of the G- function and the more generahzed H- function 
sec Mathai and Saxena (1973, 1978). 

Writing the Melhn-Barncs integral representation for the G- function and 
then simplifying the gammas by using the multiplication formula for gamma 
functions, namely, 

r(mz) = (27r)(^-'")/2m"^^-^/2r(^)r(^ + 1) . . . r(^ + ^^^), m = 1, 2 . . . 

m m 

(18) 

one has the right side of (17) simplified to the following form, 

A(-2;p,n,m) / T(—)T(l-nr + s)(pz"'/''y'ds, (19) 

n 2m JL n 

where L is a suitable contour. The integral on the right side of (19) can 
be written as an H-function, see Mathai and Saxena (1978), which can be 
reduced to a G- function due to the fact that m/n is rational. Now consider 

two independent positive real random variables x and y with the density 
functions fi{x) < for < ,t < oo, /2(y) < for < ?/ < cxd and /i(x) = 
05/2(1/) = elsewhere. Let the (s-l)st moments of x and y be denoted by 
gi{s) and g2{s) respectively. Consider u — xy. Prom statistical independence 
of X and y one has 

E{u'-') = {Ex'-'){Ey^-') = g,{s)g2{s), 

where E denotes the expected value. If the density of u is denoted by g{u) 
then from the inverse Mellin transform 

g{u) [ gi{s)g2{s)u-'ds, (20) 

where Li is a suitable contour. Now by using transformation of variables 
technique the density is given by 

g{u) = Mv)M-)v-'dv. (21) 
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But due to the uniqueness of the density of u the g{u) appearing in (20) 
and (21) must be one and the same. Let v = pt,u = pz'^^^,fi{t) = 
t-'^~"''e~*and/2(t) = e"*"''"'. Then one has 

gi{s) = r(l — nr + s), R{1 — nr + s) > 

5 =-r — >o 

n n 

and 

r h(v)h{-)v-'dv= r e-'W-'e^-'-r^dt, (22) 
Jo V Jo 

which gives 

AJz;p,n,m) =p''U—)— [ T( — )T(l-nr + s)(pz"'/'')-'ds. (23) 
n 2ni Jli n 

Comparing (23) and (19) the result is established. 

Thus in order to evaluate the basic colhson probabihty integral one has 
to represent the G-function in Theorem 3.1 in computable forms. Explicit 
computable representations of the G-function for various parameter values 
are given in Haubold and Mathai (1984a). Computable representation of a 
general G-function is available in Mathai and Saxena (1973). 

3.1 Non- resonant reaction rate 

In this case the closed-form representation of the screened nuclear reaction 
rate can be evaluated by evaluating the following integral. 

I(z;t,a,v,n,m) ^ e'^Ve'^^^^*^ ciy, (24) 
Jo 

for z > 0,t > 0,a > 0,m,n positive integers, see Haubold and Mathai 
(1986c) for details. 

The integral in (24) can be evaluated by using the following lemmas, of 
these, Lemma 3.1a can be established by going through the same process 
as in the proof of Theorem 3.1, but Lemma 3.1b needs some modifications. 
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In this case one of the random variables will be having a non-zero density 
function in the interval [0,d] and the density function will be zero outside 
this interval. Apart from this modification the derivation will remain more 
or less parallel. Hence we will list the lemmas and the theorem to follow 
without proofs. 

Lemma 3.1a. For a > 0, z > 0,n,m positive integers 

Ni(z;a,v,n,m) = / e-"Ve~^^ dy 

Jo 

Lemma 3.1b. For z > 0,d > 0,a > 0,m,n positive integers and denoting 

N2{z; d, a, v, n, m) by N2, 

N2 = [%''e-''ye-'y~"^"^dy 
Jo 

= (2.)^,„V=„-V-'Et|?LG™-;«(^|«), (26) 

n n n m 

a) = (-H \ = 

n n m 

Theorem 3.2. For a > 0,z > 0,t > 0,m,n positive integers and v a 
non-negative integer 

-N2{zi;l,ai,v - r,n,m)], 

where ai = at, z\ = zt, (^) = i)\/[r\{v — r)!], 0! — 1, iVi(.) and A^2(-) are given 
in (25) and (26) respectively. 

3.2 Resonant reaction rates 

Analytic representations of the reaction probability for the resonant reaction 
rates are discussed in detail in Haubold and Mathai (1986). The integral to 
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be evaluated has a more complicated structm^e compared to the integrals in 
Section 3.1. In this case the integral to be evaluated is where, 



_ /-oo t-exp[-at-qt-^/^] 

^'-Jo — {b - ty + ^^^^ 

for a > 0, > 0, TO, n positive integers and v a non-negative integer. This 
will be evaluated with the help of Lemmas 3.1a, 3.1b, Theorems 3.1, 3.2 and 
the following result which will be stated as a lemma. 
Lemma 3.2a. For (b - t)^ + g"^ > 0, 

/■oo - - 

[{b - tf + g\' = exp {-[{b- tf + g'\x} dx. (28) 

Now replacing the denominator on the right side of (27) by the integral in 
(28) and evaluating the double integral with the help of Lemmas 3.1a, 3.1b 
and Theorems 3.1 and 3.2 we have the following result. 

Theorem 3.3. 



= z^^z m (-1) 



k=o 9^(9^ J ki=o 

^2fe-fci^-(i;+fci+l)j^27r)(2~"-'")/2^1/2^t;+fcl+l/2 



^m+n,0 I H " I 



where Nz[.) is defined in (27). One can look at structures which are mathe- 
matically more complicated than the one in (27). For example the denomina- 
tor in (27) could be replaced by [{b — tY + g'^X' for 7 > 1. Still the techniques 
described above will work. One can replace t""/™ by (l + t)""/"^ and one can 
also look at a finite range integral going from to some number < 00 in- 
stead of to cxD. Integrals with such modifications can be tackled by using the 
combined statistical and generalized special function techniques described in 
Sections 2 and 3 of this paper. 
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